The spin-half Heisenberg antiferromagnet on two Archimedian lattices: 
From the bounce lattice to the maple-leaf lattice and beyond 



D. J. J. Farnell^, R. Darradi^, R. Schmidt-^, and J. Richter'^ 
^Division of Mathematics and Statistics, 
Faculty of Advanced Technology, 
University of Glamorgan, 
Pontypridd CF37 IDL, Wales, United Kingdom 

^Institut fiir Theoretische Physik, 
Technische Universitat Braunschweig, 
D-38106 Braunschweig, Germany 
and 

^Institut fiir Theoretische Physik, 
Otto-von-Guericke Universitat Magdeburg, 
39016 Magdeburg, Germany 



We investigate the ground state of the two-dimensional Heisenberg antiferromagnet on two 
Archimedean lattices, namely, the maple-leaf and bounce lattices as well as a generalized J- J' model 
interpolating between both systems by varying J'/ J from J' / J — (bounce limit) to J'/ J — 1 
(maple-leaf limit) and beyond. We use the coupled cluster method to high orders of approximation 
and also exact diagonalization of finite-sized lattices to discuss the ground-state magnetic long-range 
order based on data for the ground-state energy, the magnetic order parameter, the spin-spin cor- 
relation functions as well as the pitch angle between neighboring spins. Our results indicate that 
the "pure" bounce {J' /J ~ 0) and maple-leaf {J' /J = 1) Ifeisenberg antiferromagnets are mag- 
netically ordered, however, with a sublattice magnetization drastically reduced by frustration and 
quantum fluctuations. We found that magnetic long-range order is present in a wide parameter 
range < J' /J < J'c/J and that the magnetic order parameter varies only weakly with J'/ J. At 
Jc ~ 1.45J a direct first-order transition to a quantum orthogonal-dimer singlet ground state with- 
out magnetic long-range order takes place. The orthogonal-dimer state is the exact ground state in 
this large- J' regime, and so our model has similarities to the Shastry-Sutherland model. Finally, we 
use the exact diagonalization to investigate the magnetization curve. We a find a 1/3 magnetization 
plateau for J' /J > 1.07 and another one at 2/3 of saturation emerging only at large J' /J > 3. 



I. INTRODUCTION 



The study of two-dimensional (2D) quantum spin- 
half antifcrromagnetism is an interesting and challenging 
problem, in particular, if the magnetic interactions are 
frustratediii^ In 2D systems the interplay between geom- 
etry and quantum fluctuations may lead to semi-classical 
ground state (GS) phases with conventional magnetic 
long-range order (LRO) as well as to new quantum phases 
without magnetic LR04i^ The spin-i Heisenberg antifer- 
romagnet (HAF) on the eleven 2D Archimedian and re- 
lated lattices presents an excellent opportunity to inves- 
tigate the subtle balance of interactions and fluctuations 
and the role of lattice geometry.— It is well-established 
that magnetic LRO is present in the GS of the spin- 
i HAF on bipartite lattices (squaroi*^, honeycomb^"—, 
1/5-depleted square (or CaVO)^ii^, square-hexagonal- 
dodecagonal^) . However, this magnetic order can be 
weakened or even suppressed by the presence of frus- 
tration. The flrst investigation of a frustrated quantum 
HAF goes back to the early 1970s, when Anderson and 
FazekasiSiii considered the spin-i HAF on the triangu- 
lar lattice. They conjectured a magnetically disordered 
GS. However, it is now clear that there is magnetic GS 



LRO in this system, see, e.g., Refs. |4| |12M161 The pio- 
neering work of Anderson and Fazekas has formed the 
starting point for an intensive investigation of frustrated 
quantum magnetism. In particular, it has stimulated the 
search for non-magnetic quantum states in 2D magnetic 
systems. 

Another well investigated frustrated 2D system is the 
HAF on the SrCuBO lattice, which can be transformed 
by an appropriate distortion to the Shastry-Sutherland 
square lattice modeli^ with equal strength of all exchange 
bonds. Again, the GS is magnetically long-range ordered, 
see e.g. Refs. [I8I - I22I . However, the magnetic LRO may 
be destroyed by a modification of bond strengths i^^r— 

The coordination number z is quite large for the tri- 
angular (z = 6) and SrCuBO (z = 5 ) lattices, and that 
might be responsible for the semi-classical GS LRO found 
for the HAF on these frustrated Archimedian lattices. 
Thus, a "non-magnetic" quantum GS might be favored 
for lattices with lower coordination number z. Indeed, 
a regular depletion of the triangular lattice by a factor 
of 1/4 yields the Archimedian kagome lattice with co- 
ordination number z = 4. Contrary to the triangular 
lattice the GS of the spin-i HAF on the kagome lat- 
tice is most likely non-magnetic, see e.g. Refs. l 323l - [27[ 
Another frustrated model with low coordination num- 
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ber z = 3 having a non-magnetic quantum GS is the 
HAF on the Archimedian star latticei^i^"— Moreover, 
we mention that the HAF on the (non- Archimedian) 
square-kagome lattice with z = 4 has most hkely also 
a non-magnetic quantum GS4i"— The above-mentioned 
depletion of the triangular lattice by a quarter is clearly 
not the only possibility. As has been pointed out by D. 
Betts^ a regular depletion of the triangular lattice by a 
factor of 1/7 yields another translationally invariant lat- 
tice, namely the Archimedian maple-leaf latticei^i^ The 
coordination number of this lattice is z = 5 and lies be- 
tween those of the triangular (z = 6) and the kagome 
(z = 4) lattices. Moreover, there is a frustrated Archi- 
median lattice with z = 4, the so-called bounce lattice;^ 
Both the maple-leaf and the bounce lattices might be 
candidates for non-magnetic GSs. However, there are 
indications from previous studies (based on exact diag- 
onalizations of finite-sized lattices) of semi-classical GS 
magnetic LRO^i^ in the spin-i HAF on these lattices. 
This conclusion was drawn based on only two finite lat- 
tices of iV = 18 and TV = 36 sites and therefore the 
conviction in the conclusions is lessened. 

It is interesting to point out that the discussion of 
the magnetic properties of the HAF on Archimedian lat- 
tices is not only a challenging theoretical problem of 
quantum many-body physics but that it is also strongly 
relevant to experiment. Indeed, most of these lat- 
tices are found to be underlying lattice structures of 
the magnetic ions of various magnetic compounds, such 
as GaV40g^, SrCu2(B03)2^, or [Fe3(/i3-0)(MT-OAc)6- 
(H20)3][Fe3(M3-0)(M-OAc)7.5]2-7 HzO^ Recently there 
has been also synthesized the magnetic compound 
Mx[Fe(02GGH2)2NCH2P03]6 • nH20 (Ref. EO) with 
maple-leaf lattice structure as well as hybrid cobalt hy- 
droxide materials^ with maple-leaf like lattice struc- 
ture which may stimulate increasing interest in this 
lattice. Moreover, in the natural mineral spangolite 
(GueAl(SO4)(OH)i2Cl-3H20) the magnetic copper ions 
sit on the lattice sites of the maple-leaf latticei^i^ In par- 
ticular, spangolite is a very interesting magnetic system, 
since magnetic copper ions carry spin 1/2 and the ex- 
perimental data indicate that strong fluctuations at low 
temperatures are present which may prevent magnetic 
ordcring4^ In Ref. Fennell et al. propose that the 
spin-half Heisenberg model on the maple-leaf lattice with 
five different exchange integrals is the relevant model for 
this material. 

In this article we will discuss the GS properties of the 
spin-i HAF on the maple-leaf and the bounce lattices, 
see FiglH These lattices are related to each other be- 
cause the bounce lattice is equivalent to a bond-depleted 
maple-leaf lattice, see Ref. |4| and Fig. [2] Therefore we 
will consider a generalized spin-^ J -J' HAF 

H = s, • + ,/' ^ s, • ; J > ; J' > , (1) 

{ij) [y] 

where {ij) runs over all nearest-neighbor (NN) bonds of 
the bounce lattice, and [ij] runs over all additional NN 




FIG. 1: Archimedian Lattices: (a) 1/7-depleted triangular 
(maple-leaf) lattice, (b) bounce lattice. The six equivalent 
sublattices for both lattices are indicated by '1', '2', '3', '4', 
'5', '6' 



bonds present in the maple-leaf lattice, cf. Figs. [5] and [31 




FIG. 2: Illustration of the classical GS of J — J' model 
(Eq. Ul; J - solid lines; J' - dashed lines). The numbers 
1, 2, 3, 4, 5, and 6 denote the six sites in the geometrical unit 
cell. The magnetic unit cell contains three geometrical unit 
cells A, B and C. 

A well-established method that can deal effectively 
with the GS properties of infinite 2D quantum magnets 
is given by the coupled cluster method (CCM) (see, 
e.g., Refs. IH-IH and references cited therein). The 
accuracy and effectiveness of this method in relation to 
the investigation of frustrated quantum spin systems 
has been strongly improved by the implementation of 
a parallelized CCM code^ that carries out high-order 
CCM calculations. In particular, quantum phase tran- 
sitions in 2D quantum spin systems that are driven by 
frustration can be studied by this methodi ^^i^"" — 



II. THE CLASSICAL GROUND STATE 

We start with a brief illustration of the classical GS 
of the model (see Fig. ^ i.e. the are treated as 
classical vectors of length s. The classical GSs in the 
limits J' = 1 and J' = have already been discussed in 
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Ref. S Starting from the information on the classical 
GS provided there it can be found easily for the J — J' 
model. The geometrical unit unit cell (hexagon) con- 
tains 6 sites labeled by the running index n = 1,...,6. 
The magnetic unit cell is three times larger. Within 
one geometrical unit cell the pitch angle a between 
neighboring spins (e.g. spin on sites 1 and 2) around the 

hexagon is given by a = tt — arctan( ^j^, ) ■ Next-nearest 
neighbor spins on a hexagon (e.g. spins on sites 1 and 
3) are parallel. Equivalent spins in two neighboring unit 
cells are rotated by a fixed angle ±27r/3, see the spin 
directions in unit cells A, B, C shown in Fig. [2j In 
the limits J' = (bounce lattice), J' = 1 (maple-leaf 
lattice) or J' oo one has a = tt and E = —3JNs'^/2, 
a = 57r/6 and _E = -(1 + ^/2,)JNs'^/2 or a ^ 7r/3 and 
E —J'Ns'^/2, respectively. 



III. THE METHODS 

A. Coupled Cluster Method (CCM) 

We start with a brief illustration of the main features 
of the CCM. For a general overview on the CCM the in- 
terested reader is referred, e.g., to Refs. I45ll47ll60l and for 
details of the CCM computational algorithm for quantum 
spin systems (with spin quantum number s = 1/2) to 
Refs. lliliiliiiO. The starting point for a CCM calcula- 
tion is the choice of a normalized model or reference state 
1 3>) , together with a set of mutually commuting multispin 
creation operators C/ which are defined over a complete 
set of many-body configurations /. The operators C/ are 
the multispin destruction operators and arc defined to be 
the Hermitian adjoint of the C/. Wc choose {|<f>);C/} 
in such a way that we have ($|C/ = = C/|$), V/ 7^ 0. 
Note that the CCM formalism corresponds to the ther- 
modynamic limit — > 00. 

For spin systems, an appropriate choice for the CCM 
model state |$) is often a classical spin state, in which the 
most general situation is one in which each spin can point 
in an arbitrary direction. We then perform a local coor- 
dinate transformation such that all spins arc aligned in 
negative z-direction in the new coordinate frame.— i^i^ 
As a result we have 

1$) = I U; • • • ); C+ = s+, • • • , (2) 

(where the indices i,j,k,... denote arbitrary lattice 
sites) for the model state and the multispin creation op- 
erators which now consist of spin-raising operators only. 

The CCM parametcrizations of the ket and bra ground 



states are given by 

H\^) = E\^) ; (*|i? = ; 

(#| = ($|5e-^; S = l + Y.^iCY. (3) 

The correlation operators S and S contain the correlation 
coefficients Si and Si that wc must determine. Using the 
Schrodinger equation, H\'^) = E\^), we can now write 
the GS energy as = ($|e~'^i?e'^|$) and the magnetic 
order parameter (sublattice magnetization) is given by 

Ms = ^Tq^i (^kfl'^)' where sf is expressed in the 
transformed coordinate system. (Note that all magnetic 
sublatticGS carry the same sublattice magnetization.) 

To find the ket-state and bra-state correlation co- 
efficients we require that the expectation value H = 
{^\H\^) is a minimum with respect to the bra-state and 
ket-state correlation coefficients, such that the CCM ket- 
and bra-state equations are given by 

($|C7e-^i?e^|$) = (4) 
(3>|5e-^[i/, C+]e'^|$) =0 V/ 7^ 0. (5) 

The problem of determining the CCM equations now be- 
comes a pattern-matching exercise of the {CJ^ to the 
terms in e~^He^ in Eq. 

The CCM formalism is exact if we take into account all 
possible multispin configurations in the correlation oper- 
ators iS* and iS". This is, however, generally not possible for 
quantum many-body models including that studied here. 
We must therefore use the most common approximation 
scheme to truncate the expansion of S and S in the 
Eqs. ([4]) and namely the LSUBn scheme, where we 
include only n or fewer correlated spins in all configura- 
tions (or lattice animals in the language of graph theory) 
which span a range of no more than n adjacent (contigu- 
ous) lattice sites (for more details see Refs. l46H48ll50ll6Clf ) . 
For instance, one includes multispin creation operators 
of one, two, three or four spins distributed on arbitrary 
clusters of four contiguous lattice sites for the LSUB4 ap- 
proximation. The number of these fundamental configu- 
rations can be reduced exploiting lattice symmetries. In 
the CCM-LSUB8 approximation we have finally 330369 
fundamental configurations. 

Since the LSUBn approximation becomes exact for 
n — ^ C!0, it is useful to extrapolate the 'raw' LSUBn data 
in the limit n — > 00. An appropriate extrapolation rule 
for the order parameter of systems showing GS magnetic 
LRO is given by Ms(n) = oo -I- ai(l/n) -I- 02(1/71^) (see, 
e.g., Refs. Si-Iil) where the results of the LSUB2,4,6,8 
approximations are used for the extrapolation. For the 
GS energy per spin, E{n)/N = bo + 6i(l/n^) + 62(1/^'') 
is a well-tested extrapolation ansatz 4^ ~48 ,50 , 5 1 
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B. Exact diagonalization (ED) 

The Lanczos ED is a well-established many-body 
method, see, e.g., Ref. |6l|. Hence we can restrict our 
discussion of the method to some specific features rel- 
evant for our problem. The ED has been successfully 
applied to 2D frustrated quantum spin systems, see e.g. 
Refs. iHiillll. Here we follow the lines of Rcf. i. For 
the system under consideration there arc only two ap- 
propriate finite lattices, namely one with = 18 (see. 
Fig. 2 of Ref. Isil ) and another one with = 36 (see 
Fig. [3]). Note that these finite lattices do not have the 
the full lattice symmetry of the infinite lattice. Moreover 
the unit cell of these lattices is fairly large, namely it 
contains six sites. Hence, we consider the ED data as a 
complementary information to confirm or to question the 
CCM results corresponding to the thermodynamic limit 
N oo. For the finite-size order parameter we use^ 

We extrapolate the GS energy and the order parameter 
as described in Ref. (see also Refs. 16211631 and references 
therein). However, the results of these ED extrapolations 
has to be taken with caution, since they are based only 
on two data points. 



IV. RESULTS 

Henceforth, wc set J = 1 and wc consider J' as the 
active parameter in the model of Eq. ([T]). We apply high- 
order CCM up to LSUB8 and these infinite-lattice results 
are complemented by ED results of = 18 and = 36 
sites, sec Fig. [31 We choose the classical canted state 
illustrated in Sec. [IT] to be the CCM reference or model 
state. As quantum fluctuations may lead to a "quantum" 
pitch angle that is different from the classical case, we 
consider the pitch angle a in the reference state as a free 
parameter. We then determine the quantum pitch angle 
aqu by minimizing £'LSUBTi(a) with respect to a in each 
order n. 

We start with the case of the perfect Archimedian 
bounce and maple-leaf lattices, and so we set J' = 
and J' ~ 1, respectively. Results for the GS energy and 
sublatticc magnetization are given for both lattices in 
Tabic m GS energies agree well with the previously re- 
ported datai^i^ Furthermore, we confirm the previous 
findings that the GS is magnetically ordered. However, 
due to quantum fluctuations and frustration the sublat- 
ticc magnetization is drastically reduced. Using our ex- 
trapolated CCM data (see Table |l| we find that the sub- 
lattice magnetizations are only 33% of the classical value 
for the bounce lattice and 34% of the classical value for 
the maple-leaf lattice. However, these values are still 




FIG. 3: The finite maple-leaf lattice of 36 sites imposing pe- 
riodic boundary conditions used for the ED calculations. Re- 
moving the dashed J' bonds the the exchange pattern corre- 
sponds to the bounce lattice of A'' = 36 sites. 



TABLE I: CCM-results for the spin-1/2 HAFM on the bounce 
(J' — 0) and on the maple-leaf (J' — 1) lattices. E/N is the 
GS energy per spin and Afs is the sublattice magnetization. 
The LSUBn results are extrapolated using E{n)/N = feo -f 
6i(l/n^) + &2(l/n*) for the GS energy and Ms{n) = ao -f 
ai(l/n) + a2(l/n^) for the sublattice magnetization. 



bounce 


E/N 


Ms 


LSUB2 


-0.521631 


0.404343 


LSUB4 


-0.546866 


0.339357 


LSUB6 


-0.553763 


0.298249 


LSUB8 


-0.556998 


0.265252 


Extrapolated CCM 


-0.5605 


0.1657 


maple-leaf 


E/N 


Ms 


LSUB2 


-0.483470 


0.405622 


LSUB4 


-0.512309 


0.338483 


LSUB6 


-0.520378 


0.297499 


LSUB8 


-0.523861 


0.265768 


Extrapolated CCM 


-0.5279 


0.1690 



clearly above the ED estimates (which are 22% for the 
maple-leaf and 27% for the bounce lattice) of Ref. 0. We 
believe that the CCM data reported here are more re- 
liable than the ED estimates because these ED results 
were extrapolated using only two data points (A'^ = 18 
and A^ = 36)^, see also Sec. IHIBI These rather small 
values of the order parameter, which are significantly be- 
low that for the triangular lattice,— ^i^rJ^ indicate that the 
GS magnetic LRO is fragile, and one can speculate that 
slight modifications of the model parameters might lead 
to non-magnetic quantum GSs. We remark again that 
a related experimental material called spangolite^ did 
not appear to shown magnetic LRO, and that this ex- 
perimental result also spurs us on to evaluate the more 
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general model. 



z 




FIG. 4: CCM results for the GS energy per spin E/N as 
function of J'. The CCM extrapolated values in the limit n — >■ 
oo are obtained using the extrapolation scheme E{n)/N = 
ho + hi{l/nf + h2{l/nf. ED results for iV = 18 and = 36 as 
well as ED extrapolated values using the extrapolation scheme 
E{N)/N ^ xo + xiN~^^^ are also shown. The red straight 
line shows the exact energy of the orthogonal-dimer state. 
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FIG. 5: CCM results for magnetic order parameter (sub- 
lattice magnetization) A'ls as function of J'. The data for 
the quantum model are scaled by its corresponding classi- 
cal value Mf"" — 1/2. The extrapolated values in the 
limit n — i- cxj are found using the extrapolation scheme 
Ma{n) = CO + ci(l/n) + C2(l/n)^. We also show extrapolated 
finite-lattice results for the ED order parameter using the 
extrapolation scheme m^{N) = j/o + yiN~^^'^ . 



The results for the GS energy per spin are shown in 
Fig. Has a function of J'. CCM LSUBn results for the 
GS energy are clearly converging rapidly with increasing 
n for all values of J'. We have also used our ED data for 

= 18 and = 36 to extrapolate them to — oo (for 
the details of the extrapolation, sec Rcf. |^. As already 
mentioned above this ED extrapolation has to be taken 



with caution, since it is based only on two data points. 
We find a reasonable agreement between the CCM and 
ED data for the GS energy. As indicated by the ED data 
the GS energy becomes linear in J' at larger J' > 1.5. 
This is related to the existence of a dimer product eigen- 
state where all J' bonds carry a dimer singlet<^ We find 
that this singlet orthogonal-dimer eigenstate becomes the 
GS for J' > J^. Hence our model has much in common 
with the Shastry- Sutherland modelii"— that also demon- 
strates a similar exact orthogonal-dimer GS. We use the 
intersection point between the extrapolated CCM GS 
energy per site and the energy of the orthogonal-dimer 
eigenstate given by Eqd/N = —3 J'/ 8 to determine the 
transition point — 1.449. (Note that the correspond- 
ing value based on the ED data is = 1.462.) 

Next we use the CCM results for the magnetic order 
parameter Ms to discuss the stability of magnetic LRO 
as a function of J', see Fig. [SJ It is obvious that the 
magnetic LRO persists in the whole region < J' < 
J^. This conclusion is supported by the extrapolated 
ED order parameter also shown for comparison in Fig. [S] 
Interestingly the dependence of the order parameter on 
J' is fairly weak over the whole region < J' < J^. 
Thus, the extrapolated CCM order parameter varies only 
between 29% and 37% of its classical value Mf°-' = 1/2. 
This behavior might be interpreted as balanced interplay 
between increasing of frustration and increasing of the 
number of nearest neighbors when J' is growing. Our 
data for the order parameter lead to the conclusion that 
there is a direct first-order transition to the magnetically 
disordered orthogonal-dimer singlet GS. 
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FIG. 6: ED results for some selected spin-spin correlation 
functions (soSi), i — 1,2, 8, 23, for the finite lattice of A'' = 36 
sites shown in Fig. |3] The results for the quantum s = 1/2 
model are given by thin lines with symbols. The correspond- 
ing classical results are given by thick lines with the same 
color without symbols. The location of sites '0' and 'i' in 
(soSi) can be found in Fig. [S] Note that the classical curves 
for (sqSs) and (S0S23) coincide. 
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FIG. 7: Results for the quantum pitch angle function 
of J' calculated within CCM-LSUBn approximation with n — 
4 and 8. The classical result ad is shown for the sake of 
comparison. 
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FIG. 8: Results for the relative magnetization m = M/Msat 
versus applied field of strength h obtained via ED for — 36, 
where Mgat is the saturation magnetization. 



An additional confirmation of the above discussed be- 
havior comes from the ED data for the spin-spin corre- 
lation functions (siSj) presented in Fig. [51 Again we see 
that the variation of the correlation functions with J' is 
weak almost up to the transition point J^.. Moreover, for 
J' < J'^ the correlation functions of the quantum model 
behave similar to those of the classical model. 

Finally, in Fig. [7] we compare results for the classical 
pitch angle ad-, see Sec. |TT1 and for the quantum pitch 
angle calculated by the CCM . We see that both ad 
and aqu are close to each other and that there is only a 
slight variation of the pitch angle with J' (the value of 
the quantum pitch angle in CCM-LSUB8 approximation 
is aqu = TT at J' = 0, aqu ~ 0. 89577 and it is still aqu ~ 
OJUn at J' = J'J. Moreover, the LSUB4 and LSUB8 
data for aqu almost coincide. 

Altogether, our data for the order parameter, the spin- 
spin correlation functions, and the pitch angle lead to the 
conclusion that there is most likely a direct first-order 
transition from a magnetically ordered state with spin 
orientations similar to those of the classical GS to the 
magnetically disordered orthogonal-dimer singlet GS. 



V. MAGNETIZATION PROCESS 

The magnetization process of frustrated quantum mag- 
nets has attracted much attention due to the discov- 
ery of exotic parts of the magnetization curve, such as 
plateaus and jumps, see e.g. Refs. |63 - I68l The magneti- 
zation curves for the pure bounce (J' = 0) and maple-leaf 
(J' = 1) HAF were discussed already in Ref. based on 
ED data for = 36, where no indications for plateaus 
and jumps were found. On the other hand, we have al- 
ready seen that the interpolating maplc-lcaf/bounce lat- 
tice AF model considered here for larger values of J' has 




1 2 3 4 5 

J' 



FIG. 9: ED results (iV = 18, 36) for the m=l/3 plateau. Main 
panel: Width A — h2 — hi of the plateau versus J'. Inset: 
End points hi and /12 of the plateau versus J'. 



much in common with the Shastry-Sutherland model. In 
particular, that both have at zero field a orthogonal- 
dimer singlet ground state. It is well known that the 
magnetization curve of the material SrCu2 (303)2 as well 
as that of the corresponding Shastry-Sutherland model 
possesses a scries of plateaus, see e.g. Refs. I3 l69l - l73i Mo- 
tivated by this we study in this section the magnetization 
curve M{h) (where M is the total magnetization and h is 
the strength of the external magnetic field) for the inter- 
polating maple-leaf/bounce lattice AF model using ED 
for = 18 and A^ = 36 sites. ED results for the relative 
magnetization m = M/Mgat versus magnetic field h for 
A^ = 36 sites are shown in Fig. [51 In accordance with 
previous rcsults^i we do not see indications for a plateau 
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FIG. 10: ED results (iV = 18, 36) for the width of the m=2/3 
plateau in dependence on J'. 

for < J' < 1. Moreover it is obvious, that the finite- 
size singlet-triplet gap determining the size of the first 
plateau at m = is small at J' = 0, J' = 1 and J' = 1.3. 
That corresponds to our finding of a magnetically or- 
dered GS for these values of J' , and therefore the m = 
plateau should disappear for iV — > oo. However, a finite 
TO = plateau exists in that parameter region where the 
orthogonal-dimer singlet state is the zero-field GS, since 
this GS is gapped. 

Similar as for the Shastry-Sutherland model there is a 
well-pronounccd 1 /3 plateau appearing at larger values of 
J'. Interestingly, this plateau emerges already for values 
of J' below ,][. K, 1.5. In Fig. [9] we show the width of 
the 1/3 plateau versus J' (main panel) as well as the end 
points of the plateau (inset) for = 18 and N = 36. We 
observe a significant change at J' = Jp-^^^ ~ 1.07. Below 
Jp^^^ the typical finite-size behavior appears, i.e. the 
plateau width shrinks with system size N and it vanishes 
at — >■ oo. By contrast, for J' > J' there are almost 
no finite-size effects. The plateau width increases rapidly 
up to about J' w J'^ and then it remains almost constant 
for J' > J^, where the end points of the plateau grow 
linearly with J'. Moreover, from Fig. [5]wc find that for 
larger values of J' > 2 there is a jump-like transition 
between the to = and the m = 1/3 plateaus. 

As already mentioned above, for the Shastry- 
Sutherland model a series of plateaus was observed. 
For the model under consideration here we find indica- 
tions for a second plateau at to = 2/3, see Figs. [5] and 
\TU[ This plateau emerges only for quite large values of 
J' > Jp2/3 ~ ^- Again wG observe very weak finite-size 
effects of the plateau width for J' > Jp^^^ fFig. [TU)) . Fur- 
thermore, our ED data suggest an almost direct jump 
from the m — 2/3 plateau to saturation m = 1. 

Finally we have to mention, that our finite-size anal- 
ysis of the plateaus naturally could miss other plateaus 



present in infinite systems, see e.g. the discussion of the 
ED data of the m{h) curve of the Shastry-Sutherland 
model in Ref. 0. Hence, the study of the magnetization 
process needs further attention based on alternative 
methods. 



VI. CONCLUSIONS 

In this article we have treated a J-J' spin-half HAF in- 
terpolating between the HAF on the maple-leaf (J' ~ J) 
and the bounce lattice [J' = 0). Moreover, we also dis- 
cuss the GS for larger values J' > J. This antiferro- 
magnetic system is geometrically frustrated and it is re- 
lated to several magnetic materials which have been ex- 
perimentally investigated recently4Si^'^ On the classical 
level the ground state is a commensurate non-coUincar 
antiferromagnetic state. To study the quantum GS of 
the spin-half model we use the COM for infinite lattices 
and the ED for finite lattices of iV = 18 and TV = 36 
sites. 

We find evidence for a semi-classical magnetically or- 
dered commensurate non-coUinear GS in a wide range of 
the exchange ratio < J' /J < J'^/J ~ 1-5. However, 
due to frustration and quantum fluctuations the sublat- 
tice magnetization amounts about 30% of the classical 
value only. Importantly, we find that at ~ 1.45 J there 
is a (most likely) first-order transition to a magnetically 
disordered orthogonal-dimer singlet product GS which is 
the exact GS for J' / J > J'^/J. Therefore, the considered 
model is somewhat similar to the Shastry-Sutherland 
model. This similarity is also observed in the magne- 
tization curve. Based on ED data we find evidence for 
plateaus at zero magnetization for J'/J > J^/ J (which 
is related to the gapped orthogonal-dimer singlet GS), at 
1/3 of saturation for J' > 1.07 J and at 2/3 of saturation 
for J' > 3 J. The transition to the 1 /3 plateau and from 
the 2/3 plateau to saturation can be jump- like. Further 
plateaus not compatible to the system sizes iV = 18 and 
A'^ = 36 and therefore missed in the ED study may ap- 
pear in the infinite system. 

Our results may also stimulate other studies on this 
interesting 2D frustrated quantum system using other 
approximate methods to compare and contrast our 
results. In particular, the influence of modifications of 
the exchange bonds which might be relevant for real-life 
materials may lead to a destabilization of magnetic order. 
Knowing that the corresponding Shastry-Sutherland 
model exhibits a series of magnetization plateaus the 
search for further plateaus may be also an interesting 
problem to be studied by different means. 
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